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q ■ Abstract 

■ Several recent works have emphasized the role of spatial dispersion in wire media, and demonstrated that arrays 

Q ' of parallel metallic wires may behave very differently from a uniaxial local material with negative permittivity. 

Here, we investigate using local and non-local homogenization methods the effect of spatial dispersion on reflection 
' from the mushroom structure introduced by Sievenpiper. The objective of the paper is to clarify the role of spatial 

. dispersion in the mushroom structure and demonstrate that under some conditions it is suppressed. The metamaterial 

substrate, or metasurface, is modeled as a wire medium covered with an impedance surface. Surprisingly, it is found 
that in such configuration the effects of spatial dispersion may be nearly suppressed when the slab is electrically 
thin, and that the wire medium can be modeled very accurately using a local model. This result paves the way for 
the design of artificial surfaces that exploit the plasmonic-type response of the wire medium slab. 
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7j ; I. Introduction 

Artificial impedance surfaces such as the corrugated surfaces (see e.g. [1]) have been studied for decades. 
•^h ' However, after the seminal paper of D. Sievenpiper et al. [2] the interest towards the artificial or high- 
>-y impedance surfaces boomed in the literature. The exotic features and the possibility to engineer the 
^ ! response of the surface have resulted in many novel and improved applications, such as quasi-TEM [3] 
l^j 1 ; and impedance waveguides [4], [5], band-gap structures [6], low-profile antennas [7], [8], leaky-wave 
antennas [9], [10], and absorbers [1 1]— [14], just to name a few. 

In parallel with the research of new possibilities to utilize the features of these exotic surfaces, also 
0© \ attempts to fully understand the physics behind the surface behavior has been conducted. This research 
has resulted in analytical formulas for different artificial impedance surfaces providing physical insight 
and engineering tools for the designers. In some cases the analyses are based on the Floquet expansion of 
csj ■ the scattered field [15], [16], whereas some models use extraction methods to fit the model for surfaces. 
In other models the complete artificial impedance surface structure is treated as a grounded slab of a 
material with given permittivity and permeability tensors [17]. Although these models would predict the 
response of the surfaces accurately, they are either laborious to use or lack real physical basis. 

In our previous work [18] surfaces formed by an array of rectangular patches over a grounded dielectric 
slab were considered and an accurate spatially dispersive model for the patch array was derived. This model 
03 ' for the patch array was later used in the presence of a grounded dielectric slab perforated with metallic 
vias [5]. The perforated dielectric substrate was modeled in the case of electrically thin substrates as a 
uniaxial material with local negative permittivity. The corresponding artificial dielectric designated here as 
wire medium is reasonably well-known in the microwave engineering [19], [20], where the wire medium 
was initially proposed to simulate electron plasma. However, the wire medium has recently become known 
to exhibit spatial dispersion [21]. Non-local properties of the wire medium have been used favorably, for 
instance, in subwavelength imaging [22]-[24] or artificial impedance surfaces [25]-[27]. 
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In [25]-[27] the grounded wire medium slab or "Fakir's bed of nails" has been studied taking the 
spatially dispersive characteristics of the wire medium into account. These works have demonstrated that 
the wire medium may behave very differently from a uniaxial material with local negative permittivity. 
It was proven that in the limit a/h — > (h is the length of the wires and a is the lattice constant) the 
structured material behaves instead as a material with extreme anisotropy with the relative permittivity 
along the wires approaching to infinity and the relative transverse permittivity equals unity. Based on 
these results, it is natural to ask if spatial dispersion may also play an important role in the mushroom 
structure, since it can be regarded as a wire medium capped with a frequency-selective surface (FSS), as 
was done e.g. in [5], [17], [28]. In these works it was shown that the response of the mushroom-type 
artificial impedance surface may be predicted accurately by considering that the wire medium behaves 
as a uniaxial material with local negative permittivity, in apparent contradiction with the studies of [21], 
[25]-[27]. 

To study these issues, here we derive a homogenization model for the mushroom structure that fully 
takes into account the effects of spatial dispersion in the wire medium. Our objective in this paper is to 
clarify the role of spatial dispersion in the mushroom structure due to the aforementioned discrepancies 
in the treatment of the wire medium. Surprisingly, our results show that unlike in the topology for which 
the array of patches is removed (i.e. the "Fakir's bed of nails" studied in [25]), the effect of spatial 
dispersion may be nearly negligible in the mushroom structure when the wire medium slab is electrically 
thin [29]-[32]. This result is in good agreement with the results of [33], where it was also noticed that by 
connecting conductive structures, such as metal plates, to the wires the spatial dispersion in wire medium 
can be avoided. In such circumstances, the wire medium may be modeled to a good approximation as 
a uniaxial material with local negative permittivity. We demonstrate that the effects of spatial dispersion 
are significant in mushroom structures only for rather thick wire medium slabs. 

The rest of the paper is organized as follows: First, we will discuss and derive two analytical models of 
the electromagnetic response of the mushroom structure. One of the models takes into account the effect 
of spatial dispersion in the wire medium, while the other model assumes that the wire medium has a 
local response. Based on the homogenization results, we will discuss in which circumstances the spatial 
dispersion effects are suppressed. Then, in section III, we will validate the two models against full wave 
simulations, demonstrating that by taking into account spatial dispersion, we obtain an accurate model for 
the response of the mushroom structures. However, under some circumstances, that are well applicable to 
most of the practical high-impedance surfaces, the simple model neglecting the spatial dispersion in the 
wire medium is accurate and perfectly capable of predicting the response of the surface to incident plane 
waves. 



II. Analytical models for the mushroom structure 

A schematic picture of the mushroom structure comprising an array of patches over a dielectric layer 
perforated with metallic vias is shown in Fig. [Q Actually, the models proposed in this work are also 
applicable for other topologies of the mushroom structure, as will be discussed later. Here the patch array 
serves merely as a good example for its spatially dispersive electromagnetic properties have been studied 
in our previous work [18]. The surface impedance describing the electromagnetic properties of the square 
patch array reads [18]: 



with 7/ c fj = ^//i / £q£cS and the grid parameter a is given by: 



a = ^lnfshW^)), (2) 
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Fig. 1. Illustration of the mushroom structure (a) from the side and (b) from the above. The structure comprises a patch array over a 
dielectric slab perforated with metallic vias. The periodicity of the vias and the patches is a, the gap between the adjacent patches g, the 
radius of the vias ro, and the relative permittivity of the host medium Eh- 



where /c cff = koy/e c g, a is the period of the array, g is the gap between the adjacent patches, and ko is 
the wave number in free space. Furthermore, the effective relative permittivity reads: 

_ £h + 1 n\ 
£ c fi — — - — , K.J) 

where e h is the relative permittivity of the supporting host medium. 

We model the dielectric slab perforated with metallic vias as wire medium which can be described with 
an effective relative permittivity tensor having values only along its diagonal. Following the notations of 
Fig. [H the permittivity tensor for this particular case reads e v !<3 g = (u x u x + u y ii y ) + £ zz u z u z , where 
£ h is the relative permittivity of host medium and e zz is the effective relative permittivity along the vias. 
In what follows, we describe two analytical models for the electromagnetic response of the considered 
structured substrate. The main difference between the two models is related to how the wire medium 
slab is treated. The second model, described in sub-section B, considers only the effects of frequency 
dispersion in the wire medium, and assumes that the material has a local response. However, it is known 
that the wire medium suffers from spatial dispersion even at very low frequencies [21]. Also it is known 
that for an electrically thin wire medium slab (with no patch array) the charges accumulate on the tips of 
the metallic vias and the effects of spatial dispersion need to be taken into account [25]. For this reason, in 
sub-section A, we describe a non-local model that takes both TM- (with respect to z) and TEM-polarized 
waves into account in the wire medium slab. In this case we assume that the effective relative permittivity 
of the wire medium slab along the vias has both frequency and spatial dispersion. 

A. Non-local model for the wire medium 

For long wavelengths the relative effective permittivity of the wire medium along the metallic vias 
(normal direction) can be written for the non-local model as [21]: 

with k = ko^/Ft is the wave number in the host medium, q z is the z-component of the wave vector 
q = (g x , q y , q z ), and k p is the plasma wave number given as [21]: 

(k v af = — , (5) 

where a is the period of the vias (the same period with the patch array) and r is the radius of the vias 
(see Fig. [B- 

A TM-polarized incident plane wave excites TEM- and TM-polarized plane waves in the wire medium 
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slab through the capacitive array. For these two plane waves we have the following dispersion equations: 

k = ±q z , (TEM mode) (6) 

k 2 = k\ + q ■ q. (TM mode) (7) 

In [25], [26] the grounded wire medium slabs were characterized by solving the field amplitudes in all 
space. Omitting the dependency on the y-coordinate (e -jfey!/ ), we can write the field amplitudes for the 
magnetic fields similarly as: 

e jkzz _|_ p e -jkzz^ if Z > 

H x = { A TEM cos (k(z + h)) + (8) 
+Ap M cosh (7tm(^ + h)) , if —h < z < 



where k z = \Jk1 — k 2 , 7tm = y k 2 + k 2 — k 2 , and k t is the transverse wave number, which is determined 
by the angle of incidence. In order to solve for the unknowns p, Atem, and A^m, an additional boundary 
condition (ABC) is needed. In [27] an ABC was derived for the interface between air and the wire medium. 
This ABC is not applicable in the present problem, because the patch array is in galvanic connection with 
the wires. However, in [26] an ABC that properly models the connection between the wire medium and 
a metallic surface was derived, and such ABC is applicable also in our case. It should be noticed here, 
however, that this boundary condition assumes that the charge density vanishes at the connection point 
of the metallic wire and the element of the capacitive array. This condition does not always hold; for 
instance in cases where the array element size is of the same order as the diameter of the wire. Anyway, 
for high-impedance surfaces, such as the mushroom structure studied here, this is seldom the case, because 
miniaturization of the surfaces leads inevitably to large array element size with respect to the period of 
the array (and to the period of the wire medium). In our case the additional boundary condition reads on 
the wire medium side of the patch array (z = 0~) as [26]: 

^- (Lue e h u z • E + (u z x k t ) ■ H) = 
dE z , dH x 

=> k e h — k t rj —— =0. (9) 

dz dz 

As in [27], the remaining boundary conditions are obtained from the classical boundary conditions with 
the exception that in our case the transverse magnetic field is discontinuous over the patch array. So, at 
the interface between the wire medium and air (z = 0) we have the following conditions of continuous 
tangential electric and discontinuous magnetic field, respectively: 

E y \ z=Q+ - E y \ z=0 _ = 0, (10) 
#x| z=0 + - #x| z=0 - = Z~ x Ey. (11) 

Here, Z g is the surface impedance of the patch array given by (OQ). In (|9l)- (fTT|) the values for the z- 
and y-components of the electric field are calculated from ([8]) by using Maxwell's equations, taking into 
account that the material is nonlocal and that the permittivity tensor depends on the considered mode. 

Finally, using (l9l)-(fTTI) we can solve the reflection coefficient p (for the magnetic field) from ([8]) 
unambigiously: 

T1VI Tlvl 

= ^ coth M) + **f± cot (kh) + jfcZ-i - -fc 
9 £ coth (WO + ^ (kh) + ft Z~i + -fc ' 
where the relative effective permittivity along the direction of the metallic vias is written for the TM 



5 



polarization using © and © as: 

^ = -( 1 -^ F )- < 13 > 

The reflection coefficient p given in (PT21) can be used to characterize the electromagnetic properties of the 
high-impedance surface in the case of TM-polarized plane wave excitation. For a TE-polarized incident 
wave the transverse electric field does not excite the metallic vias. In this case the models presented in 
[18] can be applied for modeling of the response of the surface. Furthermore, the present model is not 
only applicable for this particular type of high-impedance surface structure but for other type of surfaces 
as well, as long as the structure consists of a wire medium slab. In these cases Z g , responsible for the 
interaction with the capacitive array, should be recalculated appropriately. 

Using (l9l)-(fT2l) we can solve also ([8]) for the amplitudes of the magnetic fields inside the high-impedance 
surface structure. We use them to determine the microscopic current on the metallic wires and use this 
quantity in order to validate our assumptions on the current phase variation in the case of the local model 
for the wire medium, as will be discussed ahead. The averaged current along the metallic wires can be 
written in terms of averaged macroscopic fields and the microscopic current, /, as (see [26]): 

<4v,z = -zl = -juje e h E z + jk t H x . (14) 



B. Local model for the wire medium 

In this model for the high-impedance surface we assume that the effective relative permittivity of the 
wire medium slab along the direction of the metallic vias can be described using the local model: 

^ = (15) 

In this approximation we have assumed the current phase variation along the wires to be minimum and the 
patch array over the wire medium slab to operate as a nearly-perfect reflector. Because of the two mirror 
boundaries the wires appear infinite for the electromagnetic wave inside the grounded wire medium slab. 
In addition, the minimum current phase variation along the vias suggests that the waves can propagate 
mainly along the transverse direction inside the wire medium (q z m 0), which makes, together with the 
infinite vias, our quasi-static approximation justified. Within this model the incident plane wave excites 
a single mode (the extraordinary wave) inside the wire medium. It should be emphasized that in general 
the properties of this extraordinary mode have nothing to do with the properties of either the TEM or 
TM waves predicted by the non-local model. This model that treats the wire medium slab using the local 
approximation for the effective permittivity will be referred to as epsilon-negative (ENG) model from here 
on in order to distinguish it from the model derived in the previous section. As discussed in the previous 
section, that model uses non-local (spatially dispersive) approximation for the effective permittivity and 
requires the use of additional boundary conditions. For this reason the model presented in the previous 
section is referred to as spatially dispersive (SD) model from here on. 

A simple but yet accurate analytical model for the mushroom structure using this local approximation 
of permittivity of the wire medium slab has been derived in our previous work [5]. In order to compare 
the results of this paper with our previous results we rewrite the results of [5] in a form similar to (fT2~l) . 
After some algebra the results of [5] can be rewritten for the magnetic field reflection coefficient in the 
following form: 

a-cothfa,ft)+%z.-- 3 fc 
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where the propagation constant along the ,2-axis in the local approximation is given as 




7ioc = \l-^-k 2 e h (17) 



When comparing the SD model with the ENG model, that is Eqs. (U2l and (1161) . we find that the 
difference between the two models lies on how they treat the surface impedance of the grounded wire 
medium slab, as stated earlier in this paper. For the SD model and the ENG model the normalized surface 
admittances read, respectively: 

Vs,sb = — coth ( lTM h) + h - cot (kh) , (18) 

7tm k 

Z/s,ENG = — COth (7lo C /i) . (19) 
7loc 

The surface admittance is given in terms of the normalized surface admittance as Y s = jy s ko/r)o. 

III. Numerical validation 

In order to compare the model derived in this paper, (PT21) . and the result of our previous work, (fT6l) , we 
use these models to characterize two particular high-impedance surfaces. We will do the characterization 
in terms of reflection phase diagrams for the electric field. We will also validate the results by numerical 
simulations obtained with Ansoft's High Frequency Structure Simulator (HFSS) [34] and with CST 
Microwave Studio [35]. Following the notations in Fig. [H the parameters for the two cases read (r = 
0.05 mm in both cases): 

1) a = 2 mm, g = 0.2 mm, h = 1mm, and Eh = 10.2, 

2) a = 1 mm, g — 0.1 mm, h = 5 mm, and = 1 

In Figs.|2a) and (b) the reflection phase according to the SD and ENG models, respectively, for the first 
example are given. Quite interestingly, the general agreement between the two models for this example 
is excellent, notwithstanding the fact that each model describes the wire medium slab with very different 
material parameters. It is also seen that both models agree very well with the HFSS results. In the absence 
of losses, the ENG model shows some spurious resonances on a very narrow frequency band in the close 
vicinity of the plasma frequency of the wire medium (/ p /y / £h = 12.1 GHz). The spurious resonances 
appear due to the fact that in the close vicinity of the plasma frequency e zz — > which according to 
(fTTI ) invalidates our assumption on an electrically thin substrate. The issue of spurious resonances has 
been discussed in more detail in [28]. The results of Fig. [2] suggest that for the considered geometry the 
effects of spatial dispersion are negligible. It should be emphasized that this characteristic is radically 
different from the results obtained for a grounded wire medium slab with no patch array, for which, as 
demonstrated in [25], the wire medium slab has a completely different electromagnetic response. 

Outside of the very narrow frequency band of the spurious resonances, both models show also an 
additional physical resonance for oblique incidence. This resonance occurs due to epsilon-near zero value 
of the e zz . 

To further understand the physical reasons why the effects of spatial dispersion are negligible for this 
example, we have plotted in Fig. @](a) the magnitude of the normalized current profile along the vias. 
It can be seen that the current varies very little along the wires. The phase of the current (not reported 
here for brevity) is also practically constant. This behavior of the electric current is consistent with the 
hypotheses used to derive the ENG model, and in particular imply that the electromagnetic fields below 
the patch grid are nearly uniform along z, i.e. d/dz « 0. Thus, in the spectral (Fourier) domain the 
electromagnetic field amplitude in the wire medium has a peak at q z « 0. This explains the suppression 
of spatial dispersion, because the nonlocal dielectric function © is coincident with the local dielectric 
function (fl~5T) when q z = 0. 
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Frequency [GHz] 

(b) 

Fig. 2. The reflection phase calculated using (a) SD model in d!2t and (b) ENG model in d!6t for the first example: a = 2 mm, g = 0.2 mm, 
h — 1mm, e r = 10.2, and ro = 0.05 mm. 

In Fig. [3] the reflection phase obtained with the SD and ENG models for the second example are 
compared against the simulation results. The models are in good agreement with each other below 20 GHz. 
The disagreement between the models becomes more noticeable as the frequency increases. The HFSS 
simulation results are in very good agreement with the results of the SD model through the whole frequency 
band. If we look at the normalized current magnitudes plotted in Fig. 0](b) for the incidence angle of 30 
degrees we see that at 10 GHz the current magnitude changes somewhat. In this case the current phase 
remains constant. However, at 25 GHz (the slab is electrically thicker) the current magnitude changes 
drastically and the current phase is no longer constant. This causes the charges to accumulate on the 
wires and the effects of spatial dispersion are no longer negligible. 

The results in Figs. [2] and [3] suggest that the SD and ENG models give the same results when the 
current along the vias is essentially constant, or equivalently when the wire medium slab is electrically 
thin (kh <C 1) and the height of the metallic pins is small as compared to the lattice constant h/a < 1. 
Indeed, it may be verified by a straightforward analysis, that under these assumptions the normalized 
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Fig. 3. The reflection phase calculated using the SD model in d 1 2b and the ENG model in d!6t for the second example: a = 1mm, 
g — 0.1mm, h = 5 mm, e r = 1, and ro = 0.05 mm. 
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Fig. 4. The current magnitudes at different frequencies (a) for the first example structure and (b) for the second example structure normalized 
to the current value on the patch array (h — 0). The reader should notice the different frequencies and scales in the plots. The angle of 
incidence is 30 degrees for all cases. 



surface admittances in (fT8b and ( fT9l ) reduce to: 

fc zz . fc zz fc h fc h nm 

y - s ^^k + ^^ = ^ (20) 

?/ s ,eng ~ -?r, (21) 

respectively. Here we have used the fact that k p ~ 1/a. 

Besides the case in which the current along the wires is nearly uniform, the effects of spatial dispersion 
may also be negligible in the regime where the wire medium is characterized by extreme anisotropy. This 
situation is illustrated in Fig. [31 where, notwithstanding the current varying drastically along the wires, 
the agreement between the ENG and the SD models is good. These results are explained by the fact 
that we operate well below the plasma frequency of the wire medium. Indeed, if we assume k p ^> k we 
have for the relative effective permittivity along the wires in the non-local model for the TM polarization 
e™ ~ and in the local model for the extraordinary mode \e l ° z c \ — ► oo =>- 7i oc = j^oV^h- Under these 
assumptions the normalized admittances in (fl"8l) and (fl9l) both reduce to: 

Vs,sd ~ y s ,ENG ~ -% cot (kh) , (22) 

which, in fact, corresponds to the normalized admittance according to the TEM model [36]. The above 
result proves that the reflection properties predicted by the two analytical models are indeed the same 
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Fig. 5. The reflection phase for (a) the second and (b) the first example for the incidence angle of 45°. The reflection phases according to 
the SD, ENG, and TEM (see [36]) approximations (plotted with solid black, dashed blue, and red curves, respectively) are compared with 
each other and against CST-simulation results (plotted with red stars). The parameters for the structures are given on page 4. 



under the regime of extreme anisotropy. 

However, when these conditions (i.e. uniform current or extreme anisotropy) are not observed, the 
response of the mushroom structures may be dominated by non-local effects. To illustrate this, we plot in 
Fig. |5] the reflection phase calculated with the SD and ENG models for the same two examples considered 
before, and an angle of incidence of 45°. The star-shaped symbols in Fig. |5] correspond to full-wave results 
obtained with the commercial simulator CST Microwave Studio [35]. As can be seen, in this example 
where h/a^> 1 the SD results concur very well with full- wave simulations, whereas the ENG model only 
gives meaningful results for frequencies well below the plasma frequency, where, as discussed before the 
material is characterized by extreme anisotropy and spatial dispersion is suppressed. Quite interestingly, 
in the regime where h/a 3> 1 the response of the structure is dominated by the effects of the TEM mode, 
consistent with the results of [24]. The TEM approximation corresponds to the case where the only mode 
excited in the wire medium is the TEM mode, i.e. the TM mode is completely discarded when calculating 
the response from the wire medium slab. According to the analysis in [36], when h/a ^> 1, the TM 
mode in the wire medium is nearly suppressed. To confirm this, we have plotted in Fig. |5] the results 



10 



obtained under such TEM approximation. These results are obtained by considering that A TM = in © 
and imposing only the boundary conditions (TTOl) — (fTTT) at the interface z — 0. The results of Fig. [5](b) (the 
plasma frequency, / p A/£h ~ 92.1GHz) show that, in fact, the TEM approximation agrees well with the 
SD model and with the simulation results whereas the ENG model begins to disagree with the simulated 
results when the electrical thickness of the slab becomes comparable with the wavelength. However, the 
results in Fig. |3a) for the first example of mushroom structure (f p / \/e^ ~ 12.1 GHz) show that in the 
vicinity of the plasma frequency the TEM model fails, as expected, whereas the SD and ENG models fit 
well with the CST simulations (until the limit ka = n, after which the homogenization fails). Indeed, the 
TEM approximation is suitable only for frequencies well below the plasma frequency. 

IV. Discussion and conclusions 

In this paper the reflection characteristics of mushroom-type high-impedance surface were studied. Two 
analytical models, namely ENG and SD models, were investigated for the electromagnetic response of 
high-impedance surfaces formed by a capacitive array over a grounded dielectric slab perforated with 
metallic pins. We derived the novel analytical (SD) model for the mushroom structure that takes the 
spatially dispersive characteristics of the grounded wire medium slab into account. Using this and a more 
simple model from our previous work (ENG model), we studied the effect of the spatial dispersion in the 
wire medium on the reflection properties of the mushroom structure. 

Surprisingly, our results demonstrate that the two models concur very well in a regime where the 
thickness of the slab is much smaller than the radiation wavelength, or alternatively when the wire medium 
is characterized by extreme anisotropy. In these conditions, the wire medium slab may, indeed, be regarded 
as a uniaxial material with negative permittivity. This property contrasts markedly with other wire medium 
topologies, for which the effects of spatial dispersion are dominant [23], [25]. This dramatically different 
behavior is explained by the fact that when the wires are connected to the metallic patches/ground plane, 
the charges no longer accumulate at the tips of the metallic vias and do not create the additional current 
flow on the vias which would cause the spatially dispersive effects. Hence, the current along the wires 
can be nearly uniform and as a consequence the polarization along the wires is out of phase with the 
macroscopic electric field, as in an electron plasma. 

Furthermore, the assumptions behind the ENG model are verified by the results of this paper. We 
validated the novel model using full-wave simulations, showing that it is very accurate even when the 
current amplitude and the phase vary considerably. We also discussed the validity limits of the ENG model. 
It can be concluded that the spatial dispersion in the wire medium of the mushroom-type high-impedance 
surface structure can be suppressed in particular designs, i.e. when the patch size is sufficiently large and 
the mushroom structure is electrically thin. This applies to most of the practical high-impedance surface 
structures. 

Hopefully the results of this paper provide new insight into the properties of mushroom structures as 
well as tools for the engineers for efficient and accurate design work. 
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